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Abstract 

We introduce in the framework of the linear approximation of Gen- 
eral relativity a natural distinction between General gauge transforma- 
tions generated by any vector field and those Special ones for which 
this vector field is a gradient. This allows to introduce geometrical 
objects that are not invariant under General gauge transformations 
but they are under Special ones. We develop then a formalism that 
strengthens the analogy of the formalisms of the electromagnetic and 
the gravitational theories in a Special relativity framework. We are 
thus able to define the energy-momentum tensor of the gravitational 
field and to fully analyze the gravitational field of isolated point masses 
or continuous distributions of them obtained by linear superpositions. 

Introduction 

This paper aims to a better understanding of the linear approximation of 
General relativity where the gravitational field is described by a space-time 
metric: 



g a /3 = + h a p (1) 

where the symmetric deviations from Minkowski's metric are assumed to be 
small, and small varying, quantities. 

From the general theory that requires that the space-time metric has to 
be a tensor under general coordinate transformations it follows that h a p are 
only denned up to a gauge transformation of the following type: 

h a f3 h af3 + d a £p + dp^ a (2) 

£ a being any vector field. 
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We define the rotational of h a p as: 

FaP/j, — dahpft — dphan (3) 

and the bi-rotational as: 

Raf3nv — —-^{df_iF a p u — d u F a p^) (4) 

or: 

Ra/3fiv — —-^paphfiv + dfs u h Q ^ — d au hp^ — dp^hav) (5) 

which is of course the linear approximation of the Riemann tensor of the 
space-time ([T]). It is well-known that only this tensor, or tensors derived 
from it, depending on second order derivatives or higher of h a p, are invariant 
under General gauge transformations. On the contrary the rotational, that 
depends on first order derivatives, gets transformed according to: 

Fapn — > F a p^ + dfj,(d a £p — dpia) (6) 

and it is invariant only in particular circumstances, as for example when the 
vector £ Q is a gradient in which case we shall say that the gauge transforma- 
tion is a Special gauge transformation. 

The main idea of this paper consists in including as geometrical objects 
of physical interest those, like F a p^ and those derived from it, that are only 
invariant under Special gauge transformations. It will be organized as follows. 
In Sect. 1 we shall write down the field equations that F a p^ satisfies following 
from its definition and Einstein's field equations in vacuum. In Sect. 2 we 
shall define the energy-momentum tensor of the gravitational field. It is 
a 2-rank tensor, quadratic in F a p^ that is conserved in vacuum and that 
unlike the 4-rank super-energy tensor has the right physical dimensions. In 
Sect. 3 we apply in detail the preceding formalism to describe in full generality 
the gravitational field of an isolated point-mass moving along a geodesic 
world-line of Minkowski's space-time. In Sect. 4 we discuss the relevance 
of several gauge fixing conditions to simplify the geodesic equations. In 
Sect. 5 we consider a rigidly rotating constant-density disk and describe the 
asymptotic behavior of the stationary gravitational field obtained by a linear 
superposition of the fields of its infinitesimal constituents. In our Concluding 
remark we consider the possibility of going a little bit beyond Einstein's 
theory. 
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1 Field equations 



We shall use the following terminology: we shall say that two potential ten- 
sors hi a j3 and h 2a i3 are Weakly equivalent if they are equivalent modulo a 
General gauge transformation. A gravitational field can be defined thus as a 
class of Weakly equivalent potentials. We shall say that two potential tensors 
are Strongly equivalent if they are equivalent modulo a Special gauge trans- 
formation. We shall call a class of Strongly equivalent potentials a First order 
implementation of the corresponding gravitational field. Any gravitational 
field can be defined either by the bi-rotational R a p\^ or of by the rotational 
F a p\ of any particular First order implementation of it. 

Let us consider any first order implementation F a p\ of a gravitational 
field. From its very definition it follows that this field satisfies the constraint: 

F[ a pp] = 0. (7) 

where [, ] is the symbol for complete anti-symmetrization, and a first group 
of differential equations, namely: 

d[ 7 F a i3]n = 0. (8) 
On the other hand using (jl]) we obtain: 

d a F^ - d„F p = -2R Plx , d a F a = R, (9) 

where Rp^ and R are the linear parts of the Ricci tensor, and the curvature 
scalar: 

R^ = R aP x^ aX R = v Pu Rpu, (10) 

and where: 

Fp = tT^m- (11) 

It follows then that: 

8 a HX = -2S^ (12) 
where Sp^ is the linear part of the Einstein tensor: 

1 

Spu = Rp^ - -Rqpn (13) 

and where: 

HaPn = F a p^ + FaTjp^ — FpT] afl (14) 
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Eqs. (JBJ) can also be written equivalently as: 

d a F a 07 = (15) 
where * is the dual operator acting on the pair of antisymmetric indices: 

F a ^= ~C 2 /l^7 ( 16 ) 
Similarly Eq. (fT2l) can be written: 

H a p]u= ^VaB-ipSu (17) 



2 The Energy-momentum tensor 

Let us introduce the 2-rank tensor: 



t a /3 — Hap/iFp Pfl + FapfiHp (18) 

that can be written also as: 

tap = 2H aptl F p fni - X -^F apX H^ x (19) 
This tensor is not symmetric in general and it is symmetric if and only if: 

F p F a/3p = 0. (20) 

We start in this section the discussion of the properties that will make of 
a tensor proportional to it: 

T a /3 = Bt a p (21) 

B being a pure numerical factor, a good candidate to be chosen as the energy- 
momentum density of the gravitational field in some particular First order 
implement at ions . 

Using the field equations (jHJ) and ffT2l we have: 

d a (H a ^F Pptl ) = l -H a ^dpF app - 2S^F Ppti (22) 

Similarly using (j!5p and (1171) we have: 

* <w * i * °w * 

d a (F H Ppp ) = -F dpH am +AS^F Pm (23) 

Using now the identity: 
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Fapfj, dp H = -F^dpHow (24) 
we obtain the following formula for the divergence of the tensor ( 1TS1) : 

d a t a p = l -(H a ^dpF app - F a ^dpH app ) + 2S^F Pptl (25) 
and remembering the definition (j!4j) we get: 

d a t«p = —dp(F a F a ) + 2F a dpF a + 2S p ^F PptM (26) 

or, simplifying: 

d a t a p = 2S^F Ppil (27) 
and therefore t a p is conserved when Einstein's field equations in vacuum: 

S a p = (28) 

are satisfied. 

Contrary to this latter property that it is independent of the First order 
implementation being used the value of the time- component: 

W(e) = t a pe a e p (29) 

where e a is any time-like vector, will depend on the implementation and has 
to be discussed in each case. We shall do it in detail in one important case 
in the next section. 

Let us consider, as a simple example, a gravitational plane-wave propagat- 
ing in the x s direction. The familiar harmonic and trace free corresponding 
potential tensor is: 

h a p = A{u){8 la 8 X p - 5 2a 5 2 p), u = t-x 3 (30) 

and the corresponding expression of the First order implementation of the 
gravitational field is: 



F ia p = \jt In A^j ((<$o 7 - S^hap - (S 0a - S 3a )h 1 p) 



(31) 



In this case we have: 



F a — 0, H ja p — F ia p, F pa pF pal3 — (32) 
and therefore the tensor (11811 is simply: 
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ta/3 — FaptjFp' (33) 

It leads to the following result: 

t a/ 3 = 2 (j^) Up, l P = d p u (34) 

and encourages the interpretation of (fTB"]) as being proportional to the energy- 
momentum tensor of the plane wave. 

3 The gravitational field of a point mass 

As early as in 1906 Poincare [I] used the Special relativity framework to 
propose a theory of the gravitational field of a point mass. In 1922 Whitehead 
[2]lj, a few years after Einstein's General relativity was successfully completed 
and tested, reconsidered a point of view similar to that of Poincare to develop 
what he thought could be a concurrent theory of General relativity avoiding 
some aspects of this theory that he disliked. 

Let us consider a point particle with mass m moving with constant 4- 
velocity u 13 along a time-like world-line W of Minkowski's space-time, and 
let x a be the coordinates of a generic event where one wants the gravitational 
field created by m to be calculated. We shall use the following notations: x a 
will be the coordinates of the intersection of the past light-cone C_ issued 
from x a with W; L a will be the null vector: 

L a = x a -x a ; (35) 

and r > will be the scalar: 

r = -u a L a (36) 

which is the space distance from x a to W. Actually we shall use the reduced 
null vector l a : 

l a = -L a (37) 
r 

instead of L a . The algebraic equalities: 

u a u a = -l, u a r = -l, i a r = o, (38) 



1 Reference [3] contains also abundant information about Whitehead's theory 
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together with the easily derived formulas below are the main side conditions 
to be taken into account in the calculations to be made in this section: 

d p u a = 0, d p r = n p , d p l a = ^(r] ap + u p l a + l p u a - l p l a ) (39) 
where d p means partial derivative with respect to x p and where: 

n p = -u p + l p (40) 

Whitehead's theory is based on the assumption that the gravitational 
field of the point-particle m is described by the following symmetric covariant 
tensor invariant by the Poincare group: 

9a0 = Va/3 —I J/3 (41) 

where we have assumed that G = c = 1. 

From Whitehead's point of view ( HTl) is the foundation of a new theory 
but in fact (jHJ) is the exact Schwarzschild solution of Einstein's equations, 
written in a particular system of coordinates that nowadays are called Kerr- 
Schild coordinates. 

This section generalizes the starting point of Whitehead while sticking 
to the point of view that we are dealing with the linear approximation of 
Einstein's theory. We assume thus Eq. (CEJ) requiring h a p to be the more 
general Poincare invariant tensor built with the geometrical quantities at our 
disposal, i.e, l a , u a and four constants A -A 3 with the dimension of length 
in the system of units that we are using: 

h a p = - (A u a Up + Airj a p + A 2 (u a lf3 + l a up) - A 3 ljp) (42) 

To maintain h a p as a Poincare invariant tensor the more general accept- 
able gauge transformation that we shall consider is that where: 

£a = -zhla + k 2 u a In r (43) 

ki and k 2 being two constants. From the third set of Eqs. (l3"9"j) it follows that 
l a is a gradient and therefore if k 2 = the transformation is a Special gauge 
transformation as defined in the Introduction. 

A General gauge transformation such as (H3l) induces the following trans- 
formation on the coefficients A -A 3 in f|42|) : 

A ^A -2k 2 , A 1 -> Ai + fci, A 2 -> A 2 + ^ + k 2 , A 3 -> A 3 + k x (44) 
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Let us consider a linear combination of the coefficients A Q -A 3 : 



C = a A + a x A x + a 2 A 2 + a 3 A 3 (45) 
Necessary and sufficient conditions that leave C invariant are the following: 

cti + a 2 + a 3 = 0, — 2a + a 2 = (46) 

Some invariant linear combinations have very special tensor meanings. For 
example we have: 

A - 2(Ai -A 2 ) = 0^R afl = (47) 

and: 

{A - 2{A l - A 2 ) = 0, and A x - A 3 = 0} & R a ^ v = (48) 
Let us consider a frame of reference such that: 

u° = 1, u* = 0, x* = 0, (49) 

Then we have: 



x 
r 

Using (j4"2"j) the components of the metric §Bj become: 



z° = i, r = -, r 2 = Y,{x k Y (so) 



g 00 = -H (51) 

r 

flto = 2 2 3 Xi (52) 

,/,, = 5ij + A ^- A ^ (53 ) 

An obvious geometrical invariant in the corresponding stationary frame 
of reference at the linear approximation is: 

Ap - A x + 2A 2 - A 3 

V-goo = 1 J" ( 54 ) 

that becomes using (|4T|) and assuming the vacuum field equations (128]) . as 
we shall do from now on unless otherwise stated: 

Ai- A 3 



V-9oo = 1 ^ — ( 55 ) 
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From it we derive the force per unit test mass at rest with respect to the 
point source: 



fi = -<W-#oo = ~ Al 2r 2 AS ( 56 ) 
leading to the forced identification: 

A x - A 3 = 2m (57) 

to guarantee the appropriate Newtonian limit. 

If Ax = Ai = (l5TI) - (l5"3l) is the Schwarzschild solution, in Kerr-Schild 
coordinates, as we mentioned before. If A 2 = A 3 and Ax = we have 
the linearized Schwarzschild solution written in curvature coordinates, and 
if A2 = A3 = and A\ = 2m we have it in harmonic, or equivalently at this 
approximation, isotropic coordinates. 

More generally it can be seen easily that any two vacuum solutions with 
the same value of m describe the same gravitational field because their po- 
tentials are necessarily Weakly equivalent. This is a simplified version of 
Birkoff's theorem at the linear approximation. 

The remaining invariant geometrical objects are the Rotation rate field, 
which is zero: 

Qij = dih 0j - djhoi = (58) 

and the linear approximation of the curvature tensor of the space metric 
above which is: 



Tfl 

Rikji = —(2(5ij8ki — SkjSii) + 3(5ijlkk + Skiklj — Salklj — Skjkk)) (59) 



Using spherical polar coordinates we may write (]5Tl)-( l53l) as the manifestly 
spherically symmetric line-element: 



_ u2 n A 2 -A z , _ , / , 2m\ , / Ax\ , M ~ 



ds 2 = - 1 )dr- 2— -drdt+ 1 + )dr z + 1 + — r'dVt 1 

\ r J r \ r J \ r J 

(60) 

Keeping in mind ( HTj) and using (142p a straightforward calculation yields: 



F a p\ = ^(2(A 2 - A 3 )l x + 2(Ax - A 2 )u x ){u a lp - u p l a ) 
~(r) Xa ({Ax - A 2 )u{3 - 2ml p) - r]\p({Ax - A 2 )u a - 2ml a )) (61) 



from which expression and (fT9"l) we obtain by a straightforward calculation: 

2K 

t a p = -^-(2ljfl - 2{l a up + hu a ) - rjap) (62) 

where: 

K = (A 1 - A 2 ){2m + A 2 - A 3 ) (63) 

As expected ( 16 2 p is not invariant under General gauge transformations. Nev- 
ertheless the gauge dependence does not come from the tensor structure, but 
from the scalar K, which actually can take any value including zero. 

Let us consider a unit time-like vector e a that we can always write as: 

e a = a u u a + a t l a + a p p a (64) 

where p a is a unit space-like vector orthogonal to u a and l a . Then we have, 
using (ED): 

T a pe a eP = ^(l + 2a 2 p ) (65) 

And in particular the time component of the tensor is, when the point mass 
is at rest: 

W(u) = T a pu a u? = ^f- (66) 
while the gravitational energy density at the Newtonian approximation is: 

It is always possible to choose a First order implementation such that W(u) = 
Wn- In fact with a Gauge transformation (l4"3l) with k 2 = A 3 — A 2 , whatever 
be the initial values, the new values of the A's are: 

A\ — > A x , A 2 ^A 3 A 3 ->A 3 (68) 
and therefore using (J37|) : 

K^(Ai- A 3 ) 2 = 4m 2 (69) 
Choosing then the numerical factor: 



we obtain the tensor: 
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T aP = ^l( 2 W/3 ~~ 2 (^ w /3 + h U a) - VaP) ( 71 ) 

that in our opinion is well qualified to be the definition of the energy- 
momentum tensor density of the gravitational field of a point mass and that, 
taken as a definition, is independent of any first order implementation. 

This tensor could have been obtained directly from the following gauge 
invariant properties: i) It is a symmetric Poincare invariant tensor built with 
l a and u a ; ii) it has the right physical dimensions; iii) it is conserved; iv) the 
corresponding energy density is definite negative and has the correct Newto- 
nian limit; v) last but not least: there is a class of first order implementations 
such that it coincides with ( [181) up to the numerical factor ( ITDi) 



4 The geodesic equations 



The geodesic equations of a test massive particle or a massless photon, in the 
gravitational field of a point mass described by the space-time metric (160]) 
are: 

d 2 x a dx^ dx v dx a 

~dX r+ ^~d\~d\ = ~d\ ^ ' 

where the function b depends on the parameter A chosen, and where the non 
zero Christoffel symbols are: 



r° - 

1 01 — 


m 










1 00 — 


m 


1 ii 


m 


n i A x - 2m 


1 33 


r 2 - 

1 12 _ 


1 

r 


At 
2r 2 ' 


r 2 - 

1 33 — 


— sin 6 cos 9 




r 3 - 

1 13 — 


1 

r 


At 
2r 2 ' 


r 3 - 

1 23 — 


cos 9 
sin 9 





(73) 

i^ 2 sin 2 fl(74) 
(75) 
(76) 

As we see, some of these symbols depend on a free parameter Ax- This 
is how the covariance of General relativity is reflected at the linear approx- 
imation. This means that full operational predictions will have to use these 
equations more than once to make it intrinsic. Most often this will require 
to use the geodesic equations above to describe the motion of test massive 
particles in conjunction with the same equations to describe the propaga- 
tion of light. In which case, in the description of the experiment, Ax can be 
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chosen at our convenience depending on each particular case. In other cir- 
cumstances, when a full operational description is impossible or too difficult 
to implement the theory may recur to some supplementary assumption. The 
simplest one when the velocity of the particles is very small is that of the 
Newtonian approximation that assumes that the gravitational field does not 
affect the propagation of light and that the coordinates x % mean what they 
mean in Euclidean geometry. But less drastic assumptions may be used. 

We see from (l5TT ) -( j53l) that if the problem suggests using Cartesian co- 
ordinates of space then A\ = 2m is very likely the most convenient choice. 
On the other hand if the problem suggests using spherical polar coordinates 
then A\ = is the choice that makes the Christoffel symbols simpler. 

5 Extended bodies: an example 

An important benefit of dealing with a linear theory and having a complete 
model of the gravitational field of a point mass is that by linear superposition 
of the gravitational field of each of its components we can gain access to a 
variety of more complex sources, without needing matching conditions nor 
prescribing the appropriate asymptotic behavior. Synge [3| considered the 
case of a homogeneous spherically symmetric source and here we shall sketch 
another example that includes rotation and non sphericity of the source. 

We consider a flat disk of radius a with constant surface mass density o 
and mass M, rigidly rotating with angular velocity uj on its fixed plane and 
center O. y % will be the Cartesian coordinates of the location Y of a generic 
element of the disk, y 3 = corresponding to the plane of the disk and y % = 
to its center; x l will be the coordinates of the location X exterior to the disk 
where the gravitational field is calculated at time t. We shall also use polar 
spherical coordinates R, 9, <p and plane polar coordinates p, a on the plane of 
the disk. The square of the space distance from Y to X will be: 

s 2 = R 2 ((sin 9 cos <p — (p/ R) cos a) 2 + (sin 9 sin — (p/ R) sin a) 2 + cos 6* 2 ) 

(77) 

The velocity components of any infinitesimal component of the disk lo- 
cated at Y will be: 

v 1 = —coy 2 = —Lupsina, v 2 = +ujy 1 = cup cos a (78) 

We shall assume now that these quantities are small of order e, say, as 
well as the quantity (p/R) 2 . Under such assumptions, neglecting the terms 
of order smaller than e, the components of u a and lp at this order are: 
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u = — 1, ui = vt, 1° — 1 + Vjri* , U = l°rii (79) 

where rii is the unit vector pointing from Y to X. At the same order we have 
also: 

J = (l + ^ cos(a -(f)) - cup sin(a - 0) ) sin 6») (80) 

With these assumptions and using the gauge fixing conditions A 2 = A 3 = 
we have, following the usual procedure to transform a discrete sum into a 
continuous distribution, the following leading terms: 



a 2 



hoo — 


2M 




R ' 






h i = 


Mua 2 n 


R 3 




ho2 — 


Mua 2 , 




h 3 = 







hij = 


2M 




R ' 





s/ , 2 (l-3cos^)j (81) 

(82) 

(83) 
(84) 

a 2 \ 

' [1 -3 cos 2 9) ) 5ij (85) 



4i? 2 



that correspond to an approximate solution with a monopole-quadrupole 
rotating source. 



6 Concluding Remark 



Maxwell equations in vacuum for a point charge e moving with arbitrary 
acceleration: 

r = £ (sa) 

where r is proper Minkowski's time, follow from the assumption that the 
electromagnetic field is the rotational of a potential of the following form: 

a = \(Biu a - B 3 (l + ri'l p )l a ) (87) 

with r and l a as defined in (refl.2) and (refl.3), £ p being the acceleration at 
the retarded point x a , and: 

Bi - B 3 = e 
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We may wonder then why we needed Einstein's equations in vacuum to 
describe the gravitational field of a point mass. The answer is that we do 
not need them, at least at the linear approximation, if we are ready to deal 
with a somewhat more general theory. In fact, if we abandon the condition 
( 1471) . from ( 1541) it follows that to maintain the Newtonian limit, we have to 
make the identification: 



A - A 1 + 2A 2 -A 3 = 2m (89) 

and then writing: 

A 1 -A 3 = 2 7 m (90) 

gamma can be considered as a second fundamental constant on which de- 
pends the gravitational field. This 7 can be identified with Edington's pa- 
rameter of the PPN formalism [I] when using the gauge fixing condition 
A 3 = 

Keep in mind though that this would lead to a difficult interpretation of 
the Einstein's tensor: 

Sap = (1 -7)— 3- (Vap + u a Uj3 - 3n a np) (91) 
where we have used the notation (|40l) . 
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